Nuclear fusion in dense matter 
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The standard theory of nuclear fusion rates in strongly interacting plasmas can be (correctly) 
derived only when the energy release, Q, is large compared to other energies in the problem. We 
exhibit a result for rates that provides a basis for calculating the finite Q corrections. Crude estimates 
indicate a significant defect in the conventional results for some regions of high density and strong 
plasma coupling. We also lay some groundwork for a path integral calculation of the new effects. 

PACS numbers: 26.20.-f 
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The calculations of nuclear fusion rates in strongly cou- 
pled equilibrium plasmas, as required for some stellar en- 
vironments, nearly all rely on one simple premise, namely 
that the rate per unit volume w s of fusion of ions Ii and 
I 2 is given by, 



= n in2 (av)K^ 2 (r = 0)/K[ 0) 2 (r = 0) , 



(1) 



where ni,na are the number densities of the species, a 
is the cross-section in vacuum and K\ t 2 is the two body 
correlator in the presence of the plasma. K[ l is the 
correlator in the absence of the plasma, and basically 
cancels out the dependence of the cross-section factor on 
the Coulomb interaction between Ii and I 2 . The calcula- 
tion of K12 has engendered a big literature, e.g. [l| and 
references cited therein, in which the correlator is cal- 
culated classically in Monte Carlo simulations at larger 
distances, leading to an effective two-body potential de- 
fined as Feff(r) = -log[X lj2 (r)]//3, where = {k B T)~ l . 
This is followed by a quantum tunneling calculation us- 
ing this potential to obtain the overlap at (near) zero 
separation. This "basically classical" approach has pro- 
vided the rates that are actually used in stellar calcu- 
lations in which the plasma is strongly coupled. It has 
been tested in a very small number of quantum path 
integral calculations |3|,[3| of the correlator under some 
rather specific conditions (very degenerate electrons, one 
component plasma). The results appear to be generally 
supportive of the basically classical approach, at least in 
some domains [J]. 

In the present paper we address the fact that the as- 
sumption (JTJ that underlies both of these approaches is 
not justified in some domains in which it is currently 
being used. A particular fusion reaction has an energy 
release Q, and it is only when Q is large that ([lj can be es- 
tablished as a valid approximation, as noted in refs.QQ. 
Here we go a step farther both in addressing the question 
of how big Q must be in order that ((TJ be usable, and 
in finding a framework for numerical evaluation in the 
cases in which it is not. The results are inconsequential 
for solar physics, but relevant in denser systems. 

For the case of two ions in and two ions out, 
Ii + 12 — > I3 + I4 we take the nuclear fusion Hamiltonian 
as a point coupling, describing the idealized case in which 



all of the energy dependence of the laboratory cross- 
section in the relevant energy range is from Coulomb 
interactions, 



H ni = ge 



drvj\ (r, t)vb\ (r, t)^ (r, t)vb 2 (r, t) + H.C, 

(2) 

where Q is the energy release in the fusion. Here the ipi 
are nonrelativistic quantum fields that describe creation 
or annihilation of the respective ions, in a Heisenberg 
picture with respect to the complete Hamiltonian. The 
fields could be Fermi or Bose; deviation from Boltzmann 
statistics for the ions is inconsequential. 

The remainder of the complete Hamiltonian is taken 
as H = Hi_2 + #3,4 + H p i + H c where H lt2 and H 3 ^ 4 
are the respective kinetic energies plus Coulomb inter- 
actions of the initial and final systems in the absence of 
the surrounding plasma, H p \ contains all of the kinetic 
energy and Coulomb interactions among themselves of 
the plasma particles, and H c is the coupling of the fusing 
particles and the fusion products to the plasma particles. 

We refer the reader to ref.Q for the derivation of the 
basic formal rate expression based on ([2]), 



dte lQt / dvU\(r,t)4(r,t)Mr,t)Mr,t) 



XV»I(0, 0)^|(0,0)^ 2 (0,0)^i (0,0) 



(3) 



where the notation indicates the thermal average 

in the medium, (O)p = Zp lr Tr[Oexp(-0H)] with Z P 
the partition function. We take a one component ionic 
plasma neutralized by degenerate electrons. In this case 
the plasma coordinates are the ionic positions Ri....R/y. 
Then we can transform the result ([3]) into, 



= C 



J dv J dt e iQt J dH 1 ,..d'R, N ,..dIi' l ,..dR! N : 



(n = r 2 = 0, Ri, ..Rtfle-P-**)* |n = r 2 = r, R'^ ..R^)i, 2 



x(r 3 = r 4 = r,R' l5 .R; 



-UH 



N 



r 3 = r 4 = 0, Ri, ..R 



N/3,4,, 



(4) 



where the first two arguments in the kets are the coordi- 
nates of the reacting ions in the fusion process and the 



2 



coordinates Ri stand for all the other ions in the plasma. 
The redundant subscripts 1,2 and 3,4 on the respective 
brackets (} are a reminder of which of the reacting ions 
are present in the states within the brackets. The multi- 
plying constant is C = Zp^nxn^g 2 ■ The individual steps 
in going from © to ((U) are simple: first the expression 
of the Heisenberg fields in terms of Schrodinger fields, 
ips{ r ) = exp[—iHt]ipH(T,t)exp[iHt], then introduction 
of the single particle states for the respective reacting 
particles, |r) = ipg(r)\0) and the explicit introduction of 
the plasma coordinates. 

The result Q is equivalent to the interaction picture 
result eq. (2.53) of ref. [|. The latter is better suited 
to perturbation expansions; the former to our present 
purposes. 

To get w s of the standard theory as given by ^ we 
strike the space and time dependence of the first bracket 
in (|4]) and take the factor involving the fusion products 
to be independent of the plasma coordinates, 

w s = C J dRi,...dR.N 

x(n =r 2 =0,R 1 ,..R A r|e-' 3H |r 1 =r 2 =0,R 1 ,..R w ) 1 , 2 

= r 4 = 0)3,4 ■ 



X f dte lQt f dr(r 3 = r 4 = r\e~ ltH3 - 4 |r 3 = r 4 = 0); 

J —OO J 



(5) 



In the expression ([5]) the first two lines give 
niri2g 2 K ly 2(0) as appear in (fTJ). The final line in (O 
calculates the phase space, Q, for the fusion products, 
and the effects of their mutual Coulomb interaction, in 
a limit in which their final energy is exactly Q and their 
total momentum exactly zero. Using g 2 VL = (av^K^}- 1 
we obtain ([T]), with the caveat that we have calculated 
the final phase space neglecting momentum and energy 
pass-through from the initial state. But this is inconse- 
quential in the large Q limit, which is for other reasons 
the domain of applicability of flSJ , as noted in ref. [E| and 
explicitly demonstrated in an example below. 

The best calculation to date of the correlator -ftTi )2 , 
implicit in the first factors of ([5]), appears to be that 
by Militzer and Pollack Using exp(— j3H) = 

[exp(— [3H/Ni)) Nl in ([5|), when Ni is sufficiently large, 
they use perturbation theory for the individual factors, 
each effectively now at high temperature, N±T. This is 
an expensive calculation because N and iVi have to be 
fairly large, and between each of the JVi factors one must 
integrate over the full manifold of r' l7 r 2 , R'^.R'^. 

To apply the same technique to the complete formula- 
tion, (dJ), we would face the additional complications of 
the time and space integrations, requiring a computation 
for each point r, t that is the equilivalent of the entire cal- 
culation of ([5|). Furthermore, as it stands, a calculation 
of (j4|) at even one point r, t is out of the question because 
of the oscillating integrands. To make progress, we look 
at an approximate form of the last bracket () in (0} for 



small values of t. For simplicity we take ions 7 2 and I4 
to be infinitely massive and to be situated at the origin; 
we can then eliminate their coordinates altogether, and 
also eliminate all mention of the center of mass position, 
r, in if?]). In the last bracket (} in (01 where the cou- 
pling is to the outgoing ions #3 and #4 we make the 
further simplification in the factor relating to the fusion 
products, 



(0|e 



— it ( H p i + H$ f 4 + H c 



10) 



A(0|e 



— UH3 4 



0) 



where 



A = 1 - itH pl - it(e 3 + e 4 )0(O) . 



(G) 



(7) 



Here A is an operator in the plasma space, with the Ri 
indices now suppressed. The label in the kets refers 
only to the position of Ii. In ([7|), <f)(r) is the operator 
for the electric potential of the plasma, which enters in 
the coupling term of the I3, I4 system to the plasma, 
He = e3(j)(rs) + e40(O). Commutators that have been 
neglected in writing ([6]) in the above form give terms of 
order t 2 and higher, and would give terms of higher order 
in Q^ 1 than those that we estimate below. 
From ([4]) the rate is now, 



w 



dt e lQt F(t)G(t) 



where 



F(t) = CTr pl 



(0|e 



-(/3-it)H 



A|0) 



(8) 



(9) 



with the trace performed in the plasma space. G is given 
by, 



G(t) 



e -i*r(*-te) 1^,(0)1 



-3/2 



(t 



-3/2 (M_\^ 
V2tt 



(10) 



In (flU|) . vj/ q '(0) is the Coulomb wave function for 
asymptotic momentum q 1 , evaluated at the origin. In 
the second line we have discarded the Coulomb interac- 
tion between ^3 and 7 4 . The function G{t) is analytic 
in the upper half t plane except for the branch point 
at ie; note that this singularity comes from the high q 2 
behavior of the integrand in (fit))) and should have the 
same form when we restore the Coulomb force between 
the final particles. A similar analysis of F(t), leaving the 
plasma out entirely, would give a structure with a sin- 
gular factor (t + i(3)~ 3 / 2 . We assume that this analytic 
structure persists in the presence of the plasma. 

Then defining the cut in the t~ x / 2 function in ([9]) to 
run from ie to ioo, and that in (t + z/3)~ 3 / 2 in F to run 
from — i(3 to —ioo, we deform the t integration contour 
to run from ioo — etoOtozoo + e and replace it with a 
simple integral of the discontinuity and a real integration 
variable, r, ending with, 



3 



w = i-^ 2 f™ dtF{t)e^ t {2TiMf/ 2 {t - ie)- 3 / 2 



= -2 



dt 



dt I 
3/2 



l « t F(i)](27rM) 3 / 2 (t-ie)- 1 /2 



1/2 (#) r*r-v^ 



F(ir) 



(11) 



the integration by parts before transforming the contour 
being required in order to avoid a non-integrable singu- 
larity at t — in the final form. 

The rate corrections of order Q _1 are now calculated 
from the linear term in the expansion of of F(t) in powers 
of t. From ^ and jTJ) wc obtain just, 



F(t) = CTtpi (0|e-^[l + itfr 1)2 ]|0> 1)2 



(12) 



the terms with t<fr(0) having cancelled through conserva- 
tion of charge, and the terms with tH p \ having cancelled 
as well. The term with unity in the final bracket gives 
the rate w s of (|5|). We estimate the contribution SF of 
the linear term in t the "basically classical" approxima- 
tion, in which plasma coordinates are eliminated in favor 
of an effective Hamiltonian H c s — KE + V e s, 



5F(ir) = n 1 n 2 g 2 C 1 T{0\e 



(13) 



where ( = f dr(r\e~^ Hett |r). We evaluate (TT3]) taking a 
minimal modification of the potential due to the plasma, 
following ref. [7j, 



SH = H eS - H 1>2 « -r/?- 1 !! ~ {r/2a) 2 ] 



(14) 



where a is the average interparticle spacing and T the 
classical plasma coupling strength. We evaluate using 



<0|ffi, 2 e- 



-/3H eH 



0} = -(0| 



d_ 

dp 



5H]e 



and 



e~ So et . 



0). (15) 



(16) 



where S = [27TT 2 (3M(Ze) 4 /4h 2 }^ 3 and £ = [T - 
(45r 3 )/(325g)], as in eq.28 of ref. 7]. Here T = Ze 2 /3/a 
where a = (3/47rn/) 1 / 3 . In the second term in (fT5|) we 
take SH at the time averaged (imaginary) time f, eq.13 
in ref.[7j. We obtain, 



5w 



dSo 

' d(3 



129 r 3 
~6AS 2 (3Q 



(17) 



The first term on the RHS of (JT7J) simply adjusts the 
phase space for the outgoing particles by adding an en- 
ergy on the order of the Gamow energy to the final state, 
note the discussion of energetics above, after (0. The re- 
maining term comes specifically from the r 2 term in the 
effective potential. (The contributions from the r inde- 
pendent modification canceled, using (d/d/3)T — T//3.) 



This r 2 term gives a fractional correction « .4MeV/Q 
for the case of the extreme conditions of 12 C + 12 C at a 
temperature of 10 8 K and a density of 10 10 gc -3 . 

A correction of the above magnitude does not necessar- 
ily create a problem for, say, the conventional calculation 
of 12 C + 12 C ^ 23 Na + p where the Q value is about 2 
MeV, and where the applications care about factors of 
ten and not about 20% corrections. But the fractional 
correction is much larger than T/Q. The next logical 
step would be to expand of the right hand () in (j4j in an 
infinite series of operators in the plasma space, of which 
(IT2l) displays the first two terms. Our conjecture is that 
the contour distortion can then be applied to give indi- 
vidual terms each of which could be calculated with the 
path integral method. 

We can get complementary information from pertur- 
bative calculations in a system in which the expansion 
in powers of Q^ 1 is not applicable. We consider a re- 
action I1+I2 — > I3, where I3 is a narrow resonance, and 
where the decay of the resonance is nearly all into chan- 
nels other than the entrance channel. The formalism 
above is applicable simply by eliminating every reference 
to ion #4 in every equation. The local interaction, ([2]), 
removing the factor, perfectly describes the limit in 
which the Brcit-Wigner formula, with Coulomb removed, 
becomes a constant times a delta function in energy. As 
context we mention the possible application to 12 C + 12 C 
fusion at temperatures in the range of a few times 10 8 K, 
where the magnitude of Q, which is negative in this case, 
is chosen to be close to the Gamow peak for the fusion 
reaction, resulting in enhancements to the rate [1]- fioj ]. 

We consider only the order in which there are two 
interactions of the distinguished ions I3 with the 

plasma; the superficial order of w in the coupling param- 
eter is e 4 but the long range part of the Coulomb force re- 
duces the order of the leading term to e 3 . For these terms 
the dimensionless strength parameter is Ai = e 2 Z 2 K]j 
with n 2 D = AimiZ 2 (3. To calculate we introduce an inter- 
action representation in which the "interaction" Hamil- 
tonian is H Cl the coupling of I\,l2, 13 to the plasma par- 
ticles, and the unperturbed Hamiltonian, Ho = H — H c . 
The calculations are fairly standard, with the replace- 
ment, e.g. within the first () in 



exp 



dt'HUt') 



(18) 



where the final subscript indicates time ordering along 
the path — i/3 to to t. The perturbation terms will 
come from the first and second terms in the expansions of 
the final exponentials of (fT8|) and its (anti-time-ordered) 
counterpart for the second () in (J4|. To distinguish the 
source of the terms that follow we introduce independent 
charges of the fusing particles and the fusion products, 
ei,e2,e3, taking e\ — e 2 = eZ = e^/2 at the end of the 
calculation. When we calculate just the terms of order 



4 



e 3 in (HJ we obtain, 

w = g^nm^C 1 JZo dt \ l ~ £ _1 «£> J{t)]e lQt J dr 
x (0, 0\ e - H ^-^ |r, ^^(rle-^IO, 0) 3 , (19) 

where 

J[t) = (ei + e 2 ) 2 (^ - f + 0t) + e 2 ^ 

-(ei + e 2 )e 3 (t 2 +i/3t), (20) 

When charge conservation is imposed, the individually 
time dependent terms in ([20]) add up to — (ei + e 2 ) 2 /3 2 /2, 
and we obtain the perturbative change in rate, 

Sw = (l/2)(ei + e 2 ) 2 KD/3w (0) , (21) 

where however the e 2 and e 2 . terms are to be dropped 
because they are exactly compensated by the changes in 
chemical potential required to keep the number densities 
of the respective species constant, as shown explicitly in 
ref [f|. In the end we obtain just the Salpeter [ll[ cor- 
rection (for Z\ = Z-2 = Z), Sw = Xiw^K Using instead 
just the term in (|20|1 in which the plasma couples only to 
the incoming ions would produce a spurious addition of 
relative order (/3Q)" 1 to our earlier large Q expansion, 
and more damaging spurious effects in a resonance case 
(where we cannot use the Q~ l expansion). 

Deriving the correct Salpeter result from a wrong equa- 
tion, ([5]), can lead one to believe that there exist further 
corrections of order Ai coming from the coupling of the 
fusion products to the plasma. Other works, e.g. refs.[9j 
12], have combined a standard screening enhancement, 
of order e r , with an resonance energy shift, of order 
r//3. In our weak coupling framework, including the lat- 
ter would be incorrect, since the fractional correction to 
order Ai is given completely by (|2ip . Thus it appears to 
be incorrect in the strongly coupled case as well, though 
the issue deserves further scrutiny. 

We have calculated the residual, non-infrared diver- 
gent, terms of order e 4 in the perturbation expansion for 
the resonance production case coming from two interac- 
tions of the plasma with the final resonance. We find 
a fractional correction of a few times A 2 L[(2/i?Gamow] 
where A 2 = e i Z 3 fr/3 5 / 2 niM~ 1 /' 2 and the function L is 
of order unity when Q is of order ficamow- Thus the 
correction is superficially of order Ti while the order Ai 
correction of (f2"Tj) is classical. However, there is also im- 
plicit h dependence in -EGamow For example in 12 C+ 12 C 
at a temperature of 10 s K , this correction would be 100% 
for densities greater than about 10 8 gc -3 . The classical 
coupling r (or A 2 ^ 3 ) has become strong for even lower 
density, so this perturbative estimate of a quantum ef- 
fect must be regarded with suspicion, and is of no value 
to phenomenology. But it is worrisome that the usual 
picture does not in any way include the physics of these 
terms, especially when applied at one hundred times the 
density at which they appear to become important. 



The most important new results in this paper are the 
presentation of the basic governing equation Q and the 
demonstration that the real time dependence therein can 
lead to big changes when Q is not sufficiently large. It 
is noteworthy that, even when we evaluated the correc- 
tion using the classical potential, the time dependence 
in the framework combines with the quantum effects in 
the short distance tunneling region to produce significant 
changes. The last, perturbative, parts of the paper serve 
as a further caveat with respect to using the existing lore 
when Q is not large. 
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